home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / EULSUM.DEM < prev    next >
Text File  |  1991-04-29  |  898b  |  38 lines
  1. PROGRAM d5r1(input,output);
  2. (* driver for routine EULSUM *)
  3. LABEL 10,99;
  4. CONST
  5.    nval=40;
  6. VAR
  7.    i,j,mval,glnterm : integer;
  8.    sum,term,x,xpower : real;
  9.    glwksp : ARRAY [1..nval] OF real;
  10.  
  11. (*$I MODFILE.PAS *)
  12. (*$I EULSUM.PAS *)
  13.  
  14. BEGIN
  15. (* evaluate ln(1+x) := x-x^2/2+x^3/3-x^4/4 ... FOR -1<x<1 *)
  16. 10:   writeln;
  17.    writeln ('How many terms in polynomial?');
  18.    writeln ('Enter n between 1 and ',nval:2,
  19.             '. (n := 0 to END)');
  20.    readln (mval);
  21.    writeln;
  22.    IF ((mval <= 0) OR (mval > nval)) THEN GOTO 99;
  23.    writeln ('x':9,'actual':14,'polynomial':14);
  24.    FOR i := -8 to 8 DO BEGIN
  25.       x := i/10.0;
  26.       sum := 0.0;
  27.       xpower := -1;
  28.       FOR j := 1 to mval DO BEGIN
  29.          xpower := -x*xpower;
  30.          term := xpower/j;
  31.          eulsum(sum,term,j)
  32.       END;
  33.       writeln (x:12:6,ln(1.0+x):12:6,sum:12:6)
  34.    END;
  35.    GOTO 10;
  36. 99:
  37. END.
  38.